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ABSTRACT 

Context. The recent discovery of a very bright type la supernova, SNLS-03D3bb (= SN 2003fg), in the Supernova Legacy Survey 
(SNLS) has raised the question of whether super-Chandrasekhar-mass white dwarf stars are needed to explain such explosions. In 
principle, such a progenitor could form by the mergers of two rather massive white dwarfs. 

Aims. Binary systems of two white dwarfs in close orbit where their total mass significantly exceeds the Chandrasekhar mass, have 
not yet been found. Therefore SNLS-03D3bb could establish the first clear case of a double-degenerate progenitor of a (peculiar) type 
la supernovae. Moreover, if this interpretation is correct, it casts some doubt on the universality of the calibration relations used to 
make SNe la distance indicators for cosmology. 

Methods. Here we present a critical discussion of the expected observational finger prints of super-Chandrasekhar-mass explosions; 
in important respects, these are not consistent with the observations of SNLS-03D3bb. 

Results. We demonstrate that the lop-sided explosion of a Chandrasekhar-mass white dwarf could provide a better explanation. 
Key words. Stars: supernovae: general - Stars: white dwarfs - Hydrodynamics - Methods: numerical 



1. Introduction 

Type la supernovae (SNe la), i.e. stellar explosions whose spec- 
tra show no hydrogen but lines of elements of intermediate-mass 
nuclei and of the iron-group, have received considerable atten- 
tion recently, mainly because of their use as distance indicators 
and as tools to determine cosmological parameters. 

Although the luminosity of nearby, well-observed SNe la 
varies by more then a factor of ten, empirical correlations be- 
tween their absolute magnitude and distance-independent prop- 
erties, such as the shape of their light curves, provide means to 
calibrate their distances. Several groups have used this approach 
to estimate the lum i nosity of SNe la at redshifts up to about 1 .5 
dRiess et al l Il998[ IPeri mutter et all 119991: iTonrv et all 120031: 
iRiess et all l2004t Ustier et al.L 12006b . Although they used dif- 
ferent calibration procedures and different distant SN samples, 
they all agree that SNe la at high redshift appear fainter than 
they should in a Universe which is matter-dominated today, sug- 
gesting that the Universe started to accelerate its expansion when 
it was about half its present age. This behavior can be interpreted 
using cosmologies that include a positive cosmological constant 
or a 'dark energy' with negative pressure. 

The empirical calibrations can be justified by the fact that 
they reduce the scatter in SN la Hubble diagrams, but an im- 
plicit assumption in all of them is that all SNe la have similar 
progenitors, namely almost spherically symmetric white dwarfs 
at or near the Chandrasekhar mass, consisting of carbon and oxy- 
gen and disrupted by a thermonuclear combustion wave. In this 
scenario the (bolometric) peak luminosity is proportional to the 
amount of newly synthesized radioactive 56 Ni ('Arnett's rule', 
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lArnetll ([1982)) and the observed peak luminosity/ lightcurve 
shape relation is the result of secondary effects such as differ- 
ent opacities and/or thermal structures. It would be difficult to 
conceive ways to obtain the observed simple relationships if the 
progenitors were to span a wide range of masses, as could be 
the case for objects formed by the merger of two white dwarfs. 
Therefore, the presently favored evolutionary track to explosion 
is accretion onto the white dwarf from a main-sequence or sub- 
giant companion star rather than a merger. 

An obvious prediction of the Chandrasekhar-mass model is 
a limiting peak bolometric luminosity of SNe la of about Mboi = 
-20, corresponding to 1 .4 M of 56 Ni, the Chan drasekhar mass . 
This is confirmed by a recent parameter study (Wo oslev et all 
120061) which shows that only extreme Chandrasekhar-mass mod- 
els can reach B-band magnitudes close to -20, i.e., models in 
which it is assumed that almost the entire star is burned. In real- 
ity, of course, the Ni mass has to be lower, allowing for at least 
about 0.3 M in intermediate-mass nuclei, as seen in the spectra 
of SNe la around maximum light, and taking into consideration 
that nuclear burning will produce stable iron-group nuclei also. 
In fact, typical deflagration Chandrasekhar-mass explosion mod- 
els produce about 0.7 M or less of radioactive Ni, and even for 
a detonation incinerating the entire white dwarf to iron-group 
nuclei, some of the ash would be in the form of stable iron and 
nickel. 

Therefore it came as a great surprise when as part of the 
SNLS on April 24, 2003, a SN la was discovered at a red- 
shift z = 0.244 with a peak luminosity in the rest frame V 
band of 20.5 mag, corresponding to an absolute magnitude of 
M v - -19.94 + 0.06 for a s tandard cosm ology dHowell et all 
120061) . Using 'Arnett's' rule (lArnetll [1982). Howell et al. infer 
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a mass of radioactive Ni of (1.29 + 0.07) M and conclude that 
SNLS-03D3bb most likely was the result of the explosion of a 
rapidly spinning super-Chandra sekhar-mass white dw arf. A sim- 
ilar conclusion was reached by Jeffer v~et~ai1 d2006l) . In fact, in 
their model the mass of the progenitor has to be about 2 M Q in 
order to account for the photospheric velocity of S NLS-03D3bb 
- abo ut 8,000 km s _1 on day +2 after B-maximum dHowell et alj 
2006; s ee the discuss i on in the supplementary material to their 
paper). iHowell et al.l (|2006) also mention briefly the possibil- 
ity of an asymmetric explosion but conclude that this effect 
would be too small to account for the high luminosity of SNLS- 
03D3bb/SN 2003fg. 

In this Letter we shall first discuss this scenario in the 
light of explosion models of super-Chandrasekhar-mass white 
dwarfs and will argue that the photometric and spectroscopic 
data seem to exclude this possibility for SNLS-03D3bb. We will 
then present an alternative scenario in which nuclear burning in 
a Chandrasekhar-mass white dwarf ignites off-center and drives 
a lop-sided explosions. Here, the high apparent luminosity of 
SNLS-03D3bb would result from a nickel blob rising inside the 
star towards the observer. Synthetic lightcurves for a toy model 
are presented which support our interpretation. In a forthcoming 
paper consequences for supernova cosmology will be discussed. 



2. Super-Chandrasekhar-mass SN la models and 
their predictions 

In principle, a close binary system of two massive white dwarfs 
can merge and form a new rapidly spinning degenerate star 
that considerably exceeds the Chandrasekhar limit in which the 
final configuration is sta bilized by rotation. It was shown by 
Miiller & Erig uchil d 1985b that such stars can have masses well 
in excess of 2 M and angular momenta sufficiently low to avoid 
the various rotational instabilities, provided that they rotate dif- 
ferentially. Additionally, it is not excluded that accretion onto 
a normal whit e dwarf might have sim ilar results under certain 
circumstances dYoon & Langerl 120051) . 

If such a star loses angular momentum or the angular mo- 
mentum is redistributed closer to rigid rotation, it will become 
unstable, contract, and ignite its nuclear fuel. Central ignition 
and off-center ignition are conceivable, as well as a spontaneous 
detonation or a deflagration. 

ISteinmetz et al.l d 19921) have carried out an extended param- 
eter study of detonations of differentially and rigidly rotating 
super-Chandrasekhar-mass C+O white dwarfs in 2D assuming 
ax i-symmetry . This i s a fair approximation since it was shown 
bv lBenz et al] d 19901) by means of SPH simulations that after a 
few revolutions the merging white dwarf pair approaches an al- 
most axially-symmetr ic hydrostatic configuration. The hope of 
ISteinmetz et alj d 19921) was that due to the lower densities of the 
rapidly spinning models, a detonation would find enough white 
dwarf material at sufficiently low densities, p < 10 7 g cirT 3 , that 
carbon and oxygen would burn to silicon and not to nickel. 

To that end, the results of the study were disappointing. The 
models did produce amounts of Ni that could easily bring them 
into agreement with very bright supernovae such as SN 199 IT 
(or SNLS-03D3bb). However, no matter whether central or off- 
center explosions were considered, or rigidly or differentially ro- 
tating progenitors of different mass, the amount of Si produced 
never exceeded a few 10~ 2 M . This is, by a considerable mar- 
gin, too little to explain the strong Si II line seen in SN la spec- 
tra, including SNLS-03D3bb. The reason is simply that the det- 
onation does not leave enough time for the star to expand and 



the amount of low-density progenitor material is small, even in 
rapidly spinning white dwarfs. Moreover, for all cases consid- 
ered the expansion velocities were much higher than those in- 
ferred for SNLS-03D3bb from the spectrum observed at day +2. 
Consequently, a prompt detonation of a super-Chandrasekhar- 
mass white dwarf is clearly ruled out. 

More recently, iPfannesI d2006l) in his thesis has com- 
puted models of both, deflagrations and detonations of super- 
Chandrasek har-mass white dwarfs . The code he used is essen- 
tially that of [ Reinecke et ail d2002|). mo dified and improved by 
iRopkd d2005l) and IS chmidt et al. (2006) for pure deflagrations 
and by iMaier & Niemeve r (2006) to handle detonation in C+O 
fuel. As far as detonations are concerned, Pfanne s's work con- 
firms the earlier results of Stein metz et alj d 1992b . Although in 
a few cases and for certain rotation laws he finds up to 15% of 
the mass is burned to intermediate-mass nuclei, barely enough to 
produce the Si II feature seen in the spectrum of SNLS-03D3bb, 
the velocity (10,000 to 20,000 km s" 1 ) of this Si would be far 
too high. 

Deflagration models of super-Chandrasekhar-mass white 
dwarfs, on the other hand, fall short of producing the amount 
of Ni required for bright SNe la. Typical differentially rotat- 
ing super-Chandrasekhar-mass models predict even less Ni than 
non-rotating or rigidly rotating Chandrasekhar-mass models. 
The reason is easy to understand. Differential rotation suppresses 
the growth of the Rayleigh-Taylor instability which is the pre- 
requisite for turbulence and, in turn, determines the rate of fuel 
consumption. Differential rotation (needed to increase the mass 
beyond the Chandrasekhar limit) thus decreases the burning rate. 
Finally, the spectra expected from such models woul d fit neither 
those of normal nor of bright SNe la dPfannesll2006h . 

The only other possible super-Chandrasekhar-mass case, 
which has not been investigated until now, is that of a 
deflagration-to-detonation transition in a rapidly spinning white 
dwarf. However, it is not clear if one would get much more 
Ni in this case than from a pure deflagration. The start of 
such a model would be the deflagration phase discussed previ- 
ously. The deflagration-to-detonation transition could occur in 
unburned C+O near the center at densities where complete nu- 
clear burning would take place. But it could also start at low den- 
sity, thus resulting in incomplete burning which would produce 
intermediate-mass nuclei and not the iron group. Which of these 
two possibilities may occur, if either, will have to be determined 
from future simulations. 



3. Lop-sided Chandrasekhar-mass explosions 

Here we suggest an alternative model to explain the luminosity 
of very bright SNe la. The model we propose has two ingredi- 
ents: first, it needs about 1 M of radioactive Ni and secondly, 
the radioactive Ni responsible for the optical display of the su- 
pernova has to be distributed unevenly on large scales. 

Numerical simulations have shown ways in which this could 
happen in a real supernova. The rather high mass of Ni could 
result from a deflagration-to-detonation transition (DDT) tak- 
ing place when the star has expanded to d ensities around 10 7 g 
cm" 3 during an initial deflagration phase jG amezo et al., 2005; 
iGolombek & Niemeveit 120051; iRopke & NiemeveA 120071) . The 
asymmetry of the Ni distribution could be an int ri nsic p rop- 
erty of a DDT at one or a few points dLivnel 1 1 999b . or 
might be a consequence of an al r eady lop-sided deflagration 
dCalder et all 12004 lR5p ke et al., |2006|) (see Fig. Q] model 
3T2d200 of R opke et al.l d2006l) as an example of a significantly 
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Fig. 1. Chemical structure of a pure-deflagration off-center 3D 
explosion of a Chandrasekhar-mass C+O white dwarf 10 s af- 
ter ignition. Half of the star is shown. Volume rendered is the 
density multiplied by the respective mass fractions: iron-group 
elements (mostly 56 Ni, shown in yellow/brown); unburned C+O 
(shown in blue). 



aspherical deflagration model; note, however, that while this de- 
flagration was sufficient to unbind the star, no DDT was as- 
sumed for this particular case). In principle, asymmetric explo- 
sions might also result if after an off-center deflagration phase 
the white dwarf stays gravitatio nally bound and a surface wave 
trigger s a secondary detona tion dPlewa et al.Ll2004tlPlewal [2006, 
but see Ropke etall 120061) . Howev er, since 'surface-detonation 
models' like those o f lPlewal ((2006) predict expansion velocities 
in excess of those infer r ed fro m the observed spectrum of SNLS- 
03D3bb dHowell et al.L H006), that particular scenario does not 
seem appropriate here. 

In order to investigate the observable consequences of a 
lop-sided distribution of radioactive nickel in a rather general 
way, we have used a class of simply-parameterised models in 
which a spherical blob of nickel-rich material is placed in a ho- 
mologously expanding spherical supernova. Our models have 
two free parameters only: the mass of Ni and its displacement. 
Guided by the numerical simulations discussed earlier, we as- 
sume that the centre-of-mass of the nickel blob is displaced from 
the origin along the +z-direction by an arbitrary velocity shift 
(Ay). Within the nickel-rich blob, a fixed initial mass fraction 
of 56 Ni is adopted (/); it is assumed that there is no 56 Ni outside 
the blob. Although simplistic, these models provide a convenient 
and readily understandable tool to quantify effects. 

In the context of SNLS-03D3bb, we will discuss one particu- 
lar realisation of the model. A slice through this particular geom- 
etry (in velocity space) is shown in Fig. [2] The model has a total 
mass of 1.4 M , an initial 56 Ni-mass of 0.89 M , a maximum 
expansion velocity of 10 4 km s _1 and uniform mass density. The 
nickel blob parameters are / = 1 and Av = 10 3 km s . Wider 
discussion of these lop-sided models and their implications will 
be presented by Sim et al. (in preparation). 

U sing th e Monte Carlo radiative transfer code described by 
ISiml d2007l) . bolometric lightcurves were computed for various 
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Fig. 2. Slice through the y-z plane showing the distribution of 
56 Ni adopted in the model. The dark grey area indicates the re- 
gions in which there is no 56 Ni while the light grey area, the 
centre of which is displaced by Av = 10 3 km s from the ori- 
gin, contains pure 56 Ni (/ = 1). The models is symmetric under 
rotation about the z-axis. 
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Fig. 3. Bolometric light curves computed from the model. The 
dashed line is the angle-averaged light curve while the solid his- 
togram indicates the light curve as viewed from the +z-direction. 



viewing angles to the model. A uniform, grey-absorption cross- 
section of 0.1 cm 2 g -1 was adopted; this simple treatment is ex- 
pected to maximise the apparent viewing-angle dependence of 
the light curve and so allows a limit on the scale of any such ef- 
fect to be established. In Fig. [3] two of the computed light curves 
are shown: first, the angle-averaged light curve (dashed line; this 
is the lightcurve obtained by averaging over those seen from a 
large number of randomly oriented lines-of-sight) and secondly, 
the light curve recorded by a distant observer on the -i-z-axis (his- 
togram). 

As one would expect, the angle-averaged light curve from 
the model cannot account for the observed peak brightness of 

SNLS-03D3bb; the peak magnitude is only 19. 58 mag com- 

pared to M\, Q \ = -19.87 + 0.06 mag as reported bv lHowell et alj 
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(2006). However, the lop-sided nature of the model causes the 
light curve to be dependent on viewing angle such that if viewed 
along the direction in which the nickel blob is displaced, the peak 
is nearly 0.3 mag brighter. This effect brings us to consistency 
with the observed peak brightness of SNLS-03D3bb. 

A valid question, of course, is whether our model meets the 
velocity constraint of SNLS-03D3bb. Of course, we cannot an- 
swer this question for a "real" 3D off-center DDT of the sort we 
propose because hydrodynamical simulations for such a model 
have not yet been performed. Most likely the velocity will be 
similar to what is obtained in publ ished (less asymmetr ic) DDT 
models. For example, model (a) of Ga mezo et al.l d2005l) predicts 
a maximum expansion velocity of the outermost layers of about 
10,000 km s -1 for a total Ni mass of 0.8 AT©. On day +2 (when 
the spectrum of SNLS-03D3bb was taken) and later, most SNe 
la show constant Si velocity, indicating that the Si-layer is out- 
side the photosphere. For the Gamezo et el. model, this means 
that the Si velocity at that epoch would be ~8,000 km s , in fair 
agreement with SNLS-03D3bb. 

We can be a little more precise for the toy model we used 
here. We assumed a maximum velocity of 10,000 km s _1 (this 
is somewhat arbitrary, however, and in principle we could have 
adopted a lower or a higher velocity). As in the DDT model 
discussed above, such a peak ejection velocity does not cause 
a problem because the majority of the ejecta, even in the nar- 
row region above the Ni-bubble on the front hemisphere of the 
model, has lower velocities than the maximum value. Also, pro- 
jection effects will always cause the line-of-sight velocity de- 
duced from the spectrum to be smaller than the true velocity. In 
fact, most of the material in front of the Ni-bubble has a typi- 
cal line-of-sight velocity of 6,000 km s _1 only. Thus, although 
we cannot make detailed predictions for the spectrum on day +2 
without undertaking considerably more sophisticated radiative 
transfer calculations , we conclude that the velocities reported by 
iHowell et al. (2006) are not inconsistent with our model. 

The model considered here is tuned in a number of re- 
spects. First, the combination of Ni mass (0.89 M Q ) and geom- 
etry (Fig. |2]i is arranged to produce close to the maximum ef- 
fect possible. Further increasing the mass of the Ni blob does 

not allow a peak magnitude brighter than 20 to be obtained 

from any direction; although more Ni makes the angle-averaged 
lightcurve brighter, it does this at the expense of suppressing 
the angular variation since the mass of the surrounding shell be- 
comes smaller. For lower Ni masses it is possible to obtain more 
significant angular variations but the absolute magnitude rapidly 
becomes too faint. 

Secondly, the treatment of opacity is also extreme. The 
adopted cross-section (0.1 cm 2 g" 1 ) may be optimistically high 
for material in which significant nuclear burning has not taken 
place. Since it is the opacity in the outer regions which leads to 
the angular variations in the lightcurve, it is likely that even these 
outer regions much be burned to at least intermediate mass ele- 
ments. Conversely, the adopted cross-section may be too low for 
the region inside the Ni-rich blob; if the opacity were higher, the 
lightcurve would rise more slowly and have a correspondingly 
dimmer peak. To investigate the effects of the opacity treatment 
in detail would require much more complex non-grey, non-LTE 
radiation transport which lies beyond the scope of this letter but 
which is vital for fully understanding these effects in the future. 

4. Conclusions 

In light of these considerations, we draw the following conclu- 
sions. It may be possible to observe lightcurves with peak bright- 



ness close to that of SNLS-03D3bb/SN 2003fg from supernovae 
with a lop-sided distribution of radioactive material, and the to- 
tal mass of such material is only > 0.9 M Q ; thus this model is a 
viable alternative to the proposed super-Chandrasekhar theory. 

However, since the conditions on a simple off-centre explo- 
sion model which can explain the brightness of SNLS-03D3bb 
are quite tight, we can make several predictions from this hy- 
pothesis which may be testable via future observations. First, 
such objects must be rare since both a moderately high nickel 
mass (> 0.9 M Q ) and rather special viewing direction are re- 
quired. Secondly, SNLS-03D3bb lies very close to the limit of 
brightness that can be readily achieved - thus any supernovae 
more than one or two tenths of a magnitude brighter could not 
be explained in this way. Thirdly, the asymmetries postulated in 
our model should leave markers in the nebular spectra, such as 
a net redshift of the lines of intermediate-mass elements. And 
finally, our calculations suggest that extremely bright SNe la 
should have shorter-than-average light-curve rise times. 
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